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We study, both experimentally and theoretically, short-time modifications of the decay of excita- 
tions in a Bose-Einstein Condensate (BEC) embedded in an optical lattice. Strong enhancement of 
the decay is observed compared to the Golden Rule results. This enhancement of decay increases 
with the lattice depth. It indicates that the description of decay modifications of few-body quantum 
systems also holds for decay of many-body excitations of a BEC. 



Relaxation and decoherence of excited quantum states, 
caused by their coupling to the environment, is a ubiqui- 
tous phenomenon. Yet it still poses a fundamental ques- 
tion: can we describe its effects on complex, many-body 
quantum states, as in simple quantum systems? In few- 
body systems, control by fast modulation or short-time 
duration may result in either relaxation slowdown (quan- 
tum Zeno effect - QZE) or its speedup (anti-Zeno effect - 
AZE) [11,1,13 compared to the Golden Rule (GR) re- 
sult [5[. Here we experimentally and theoretically study 
these concepts for the decay of excitations in a Bose- 
Einstein Condensate (BEC). To this end, we apply the 
control methods of short-time duration and reservoir en- 
gineering [f| 0] to momentum excitations of a BEC in a 
deep one-dimensional (ID) optical lattice. 

We first use two-photon Bragg spectroscopy [1, 1] to 
characterize the excitation spectrum at different lattice 
depths. We then measure the decay rate of these exci- 
tations, a process describable as the Beliaev damping of 
Bloch-Bogoliubov states [l(| [H| . We find the decay rate 
to increase with lattice depth. This increase is in con- 
tradiction with the long-time GR theory [ijj]. We show 
that the non-exponential (AZE) regime predicted by the 
universal theory of short-time decay , and hitherto ob- 
served for single atoms in optical lattices Q, is in good 
agreement with the experimental results for this many- 
body process. 

Our experiment involves a BEC of - 300, 000 87 Rb 
atoms in state \F — 2,mp = 2), held in a cigar-shaped 
magnetic trap, with trapping frequencies uj±_ = 350 Hz, 
uj z = 30 Hz. The chemical potential is ~ 4 kHz, and 
the temperature is smaller than 100 nK. To alter or en- 
gineer the excitation spectrum of the BEC, we adiabat- 
ically load it onto a ID optical lattice [13], created by 
two far-detuned (0.5 nm) laser beams with wavevector 
fei, counter-propagating along the z direction (Fig. [Ha)). 
Our adiabatic ramp-up time of 250/xs is much longer than 
the timescale of the adiabaticity condition dVi a t/dt << 
l6E\/h [3, which is - 125/iS for V lat ^ 50E R . This 
ensures adiabaticity with respect to higher band excita- 
tions, while the many-body dynamics and trap dynamics 
are nearly frozen-out. These conditions allow us to study 
the dynamics of momentum excitations of a condensate 
adiabatically loaded onto a ID lattice. In the Supporting 
Information we quantitatively analyze the non-adiabatic 
effects as a function of lattice depth, following which we 
limit our experiment to lattices of Vi a t < 50-Er. Once 



FIG. 1: (a) A momentum-space absorption image (38 ms 
time-of-fiight) of the adiabatic switching on and off of an op- 
tical lattice of depth 5Er. Similar images are obtained for 
deep lattices as well (Vi a t — 50Er). (b) A 38 ms time-of- 
flight image of the Bragg excited BEC with k — 0.9fci, inside 
a lattice of 5Er. 



the BEC is thus loaded into the lattice ground state, we 
coherently excite it close to the Brillouin-zone edge by 
a pair of Bragg beams (detuned from resonance by 0.2 
nm) . This two-photon Bragg spectroscopy d, 1] allows 
us to excite and probe the BEC at a given momentum 
k = 0.9fci (along the z direction) and energy The 
time-of-flight image after the Bragg excitation (Figure 
QJb)) shows the condensate, the excitation and the de- 
cay products (occupied modes with \q\ < |fc|), created 
through collisions between the BEC and the excitation. 

The perturbative excitations of the BEC are collec- 
tive (many-body) quasi-particles in the optical-lattice 
Bloch-state basis. They can be analyzed by 3D Bloch- 
Bogoliubov equations [101 ] in terms of the condensate 
wavefunction ^ , embedded in a periodic potential in the 
axial direction, and assuming planar excitations in the 
transverse (_L = x,y) directions 18]: 
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Here m is the mass of a 87 Rb atom, Vi a t is the depth of the 
optical lattice, g is the interaction strength, n is the 3D 
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density of the BEC, k± is the transverse momentum and 
Uj q and Vj q are the amplitudes of the atomic components 
of a Bogoliubov excitation in Bloch band j with quasi- 
momentum q. Solving Eqs. ([1]) self-consistently, we find 
the excitation spectrum ujj(q,k±), and the amplitudes 
Ujq(z) and Vj q (z). 
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FIG. 2: Excitation lineshapes measured as a function of fre- 
quency using a 2ms Bragg pulse of k = 0.9fci, without a 
lattice (a) and in an optical lattice of depth Viat = WEr (b). 
The averaged data points are shown in blue, with error-bars 
depicting the error of the mean. The red lines are gaussian 
fits to the experimental data. 

In Fig. [2] we plot the measured excitation lineshape 
for k = 0.9fci as a function of the excitation frequency lo, 
with and without an optical lattice potential. Without 
the lattice, the lineshape has a single peak at the fre- 
quency of the Bogoliubov excitation (Fig. HJa)). In the 
presence of the lattice, it exhibits a double-peak struc- 
ture, centered at the frequencies of the Bloch-Bogoliubov 
excitations of the first and second bands, as calculated 
from Eq. ([I]) (Fig. OJb)). For deep lattices the bands 
are spectrally separated, and can therefore be specifi- 
cally addressed using Bragg spectroscopy. We note that 
these lineshapes indicate a qualitative difference in the 
response of the system as a function of lattice depth. 
This does not affect the results of this paper due to the 
normalization method used (see below). 

In Fig. [3]we plot the calculated Bloch-Bogoliubov exci- 
tation spectrum for lattices of different depths, and trans- 
verse momentum k±_ = 0. For weak lattices of depth 
Viat ~ Er (Er = Ti 2 k 2 L /2m is Er the recoil energy), the 
lowest band spectrum resembles the standard Bogoliubov 
spectrum: a linear phonon branch at low momenta that 
becomes parabolic at high momenta 0, [l4| . As the lat- 
tice becomes deeper, for Vi at > Er, the energy splitting 
near the Brillouin zone edge curves the parabolic spec- 
trum toward a constant value, eventually making it com- 
pletely convex at Vi a t > ^Er. Correspondingly, the width 
of the lowest band decreases exponentially, according to 

the relation [l(| huji(k L ) cx E R e~ 2 ^ Vl ° t/ER . 

The significant narrowing of the first band as a function 
of lattice depth is clearly visible in Fig. Ojnset), in which 
we plot the calculated and measured excitation energy 
for k — 0.9fci in the first band as a function of lattice 
depth. In the following, we concentrate on excitations in 
the first band, and therefore only study lattices deeper 
than Viat = %Er. This allows us to spectrally separate 
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FIG. 3: Calculated Bloch-Bogoliubov spectrum of the first 
and second bands for different lattice depths (Vi a t = OEr - 
blue dashed-dotted, Viat = IEr - red dotted, Vi a t = 15Er 
- magenta solid). As the lattice depth increases, the energy 
band flattens out exponentially. Inset: measured (red) and 
calculated (green) resonant energy of the first-band k — 0.9/cl 
quasi-momentum excitation as a function of lattice depth 
(compare to the vertical line in the main figure). 



the bands using Bragg spectroscopy. 

The coherent quasi-particle mode created by Bragg 
spectroscopy interacts with the BEC and decays into the 
continuum of other quasi-particle modes. The decay of 
this many-body coherent excitation into pairs of quasi- 
particle modes forming a continuum is governed by the 
Beliaev three- wave-mixing Hamiltonian [9( , and is known 
as Beliaev decay [IJ [ll| . 

The decay of the excitation is measured using the pulse 
sequence depicted in the inset of Fig. [4ja) . The lattice 
beams are turned on and off adiabatically over 250/is, 
and kept fixed for 500/iS. A Bragg excitation pulse at a 
frequency difference of 6kHz and Rabi frequency of 1kHz 
is applied during the first lOOus after the lattice turn-on, 
allowing for 400/xs of decay As a normalization, the 
experiment is repeated with a Bragg pulse applied dur- 
ing the last 100/j.s of the lattice existence, such that no 
time is left for decay. The ratio of the number of atoms 
remaining excited in the first experiment (Pf) to that 
recorded in the second one (Pi), measured from time-of- 
flight images such as Fig. [T]d [15| (taken by switching the 
trap off immediately following the pulse sequence and al- 
lowing 38 msec of expansion) gives the surviving fraction 
of coherently excited atoms f c (t)- 



fc(t) = 



P f (t) _ P L e 



Pi 



P, 



- (2) 



Here R(t) is the time-dependent rate of decay. This nor- 
malization method rejects "common-mode" noise in the 
experiment. 

Figure IHa) presents the measured fraction of coher- 
ently excited atoms that did not decay, as a function of 
lattice depth. It can be seen that the decay rate increases 
with the lattice depth. As shown below, this result com- 
pletely violates the Golden- Rule (GR) behavior [l2|, but 
agrees with a decay formalism that includes the effects 
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of the short-time non-exponential regime. 

The decay rate R of a state whose unperturbed energy 
is E into a reservoir is usually calculated using the Fermi 
Golden Rule 0, 



2tt 

Rgr (E ) = —p(E )\V(E )\ 2 , 



(3) 



where p(E) is the reservoir density of states and V(E) is 
the system-reservoir coupling matrix element. The GR 
presumes the Markov limit of observation times much 
longer than the correlation (memory) time t c of the reser- 
voir response. 

The change in the excitation spectrum as a function of 
lattice depth affects the decay rate in a manner that can 
be understood in terms of momentum and energy conser- 
vation of the decay process. The excited q = 0.9/cl mode 
can decay into a pair of modes, such that the sum of their 
momenta and energy equals the momentum and energy 
of the excitation. A concave spectrum, as in weak lattices 
{Viat 5 4.Er, see Fig. [3]), implies that the decay products 
have less energy than the excitation, and therefore decay 
is possible, as excess energy is carried by transverse ex- 
citations. By contrast, in deep lattices (Vi a t > 4Er), the 
completely convex spectrum (Fig. [3]) implies that the ex- 
citation has less energy than any pair of decay products 
12]. In this scenario the excitation can be described as 
a quasi-particle with dissipationless flow, analogous with 
the maxon quasi-particle of liquid He [16j |. This effect 
completely inhibits the coupling between the excitation 
and the available modes for decay at long times (compat- 
ible with the GR - see black dashed-dotted line in Fig. 

However, for finite-time decay, for which the energy- 
conservation constraint is relaxed through time-energy 
uncertainty, the above GR-based results are no longer 
valid. Since the lowest band narrows down exponentially 
with lattice depth (Fig. [3]), for any finite duration of an 
experiment there is a lattice depth beyond which the de- 
cay rate calculated by the GR (Eq. $5§) strongly deviates 
from the actual result. In our system the reservoir corre- 
lation time t c is exceptionally long (hundreds of microsec- 
onds), so that the integrated decay probability over the 
experimental time duration drastically deviates from the 
GR exponential decay [l7|. The timescale of the reser- 
voir correlation time is dictated by the bandwidth of the 
available modes for decay. In our system this bandwidth 
is given by the width of the lowest Bloch band, which is 
on the order of a few kHz. 

The above account can be given in terms of a more gen- 
eral ("universal") expression than Eq. © 2]. The uni- 
versal expression allows for arbitrary observation times, 
including t << t Cl and hence the entire spectral depen- 
dence of the excitation and the reservoir response: 



/>oo 

R Eo {t) = 2tt F{lu~ E /fi, t) G(u)du 
Jo 



(4) 



volution) of two spectral functions: F (uj — Eo/h,t), the 
time-dependent spectral function of the excitation cen- 
tered at the unperturbed energy, Eq, and G(u>), the re- 
sponse spectrum of the reservoir, which is the Fourier 
transform of its non-Markov memory functional &(t) that 
dies out at t > t c . Control of the decay is achievable by 
altering this overlap. Here we explore this simple expres- 
sion for an interacting many-body system, as compared 
to the non- interacting atoms studied before 3 . We com- 
bine reservoir engineering of G(oj) with system probing 
at short times t < t c that affects F (u> — E/h), so as to 
achieve significant changes in the overlap of the two func- 
tions and therefore of the decay. 

The reservoir spectral response function can be written 



G{uik,q) = p{uk,q)\Ak,q\' 



(5) 



where j4fc,q are the Beliaev matrix elements given in 
terms of the Bogoliubov amplitudes Uj q and Vj q from Eq. 

is the density 

of states available for decay calculated from Eq. (fTJ) • 
As the lattice deepens, the Beliaev matrix elements in- 
crease, and the spectral response of the reservoir is sig- 
nificantly enhanced. However, in the infinite-time limit 
associated with the GR (Eq. ([3|)), the system's spectral 
excitation function is 



F(E k /h - uj k>q , t — > oo) = 5(E k /h - u>k,, 



(6) 



This expression relates the decay rate to the overlap (con- 



i.e. it becomes localized at E k /Ti, the resonant frequency 
of the excitation. As the lattice depth increases, this ex- 
citation peak is pushed to lower energies, such that the 
spectral overlap between Eqs. (O and ^ goes to zero. 
Therefore, at the GR limit the decay rate of the exci- 
tation is expected to be completely suppressed above a 
certain lattice depth By contrast, in our experiment, 
the spectral excitation function of the system is spread 
around E k /fi 

F(E k /h - u k , q , At) - {Eh/n _ Uhtq)2At , (?) 

with a width (1/At) given by the inverse duration of the 
experiment. The GR result fails if this width becomes 
comparable to the energy spread of the modes available 
for decay 0], as in our experiment. 

In Fig. Ela) we show the theoretical prediction 
(dashed-red line) for the surviving coherently-excited 
fraction of atoms f c , as a function of lattice depth, along- 
side the experimental results. This is based on numerical 
integration of Eq. $2$ using Eqs. (j4|), ([5|) and (|7|), over 
the experimental time interval (400/is). The theoretical 
curve agrees with the experimental results and stands in 
contrast to the GR limit, calculated by integrating Eq. 
© over 400/Lts with R(t) = R GR (dashed-dotted black), 
which predicts that there is strictly no decay for the lat- 
tice depths of Vi a t ~ 4Er in Fig. [U(a). Both theory and 
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FIG. 4: (a) Measured and calculated surviving fraction of co- 
herent atoms, as a function of lattice depth. The averaged 
experimental data points and errorbars (depicting the error 
of the mean) are plotted in blue. The solid blue line con- 
necting the experimental points serves as a guide to the eye. 
The data clearly shows decay enhancement for deep lattices. 
The calculation is given by the dashed-red curve (see text for 
details). The dashed-dotted black curve depicts the GR re- 
sult, showing no decay for Viat > 4_Er. Inset: Sequence of 
excitation and decay measurements which gives the surviving 
fraction of coherent atoms f c (t) (see text). Following this se- 
quence the magnetic trap is switched off at t — 1000/isec, and 
a time-of-flight image is taken following 38msec of free expan- 
sion, (b) Calculated surviving fraction of coherent atoms as 
a function of time, for 3 different lattice depths (Vi a t = 5Er 
- solid blue, Viat = 10 Er - dashed red, Viat = 20Er - dashed- 
dotted black). 



experiment indicate that, as the lowest band of the spec- 
trum flattens out with increasing lattice depth, finite- 
time effects dominate to cause substantial enhancement 
of decay, as compared to the infinite-time (GR) decay. 
This is the central result of this paper, indicating that 
Eq. ((D) applies to our experiment. 

In Fig. we plot the calculated time dependence of 
the decay for different lattice depths. For all of these lat- 
tice depths GR theory predicts strictly no decay. These 



results show that the main part of the decay occurs in the 
first 100/is, which are difficult to access experimentally. 
We have also measured the remaining coherent fraction 
as a function of lattice depth for decay over 100/iS and 
250/xs, and have found results similar to Fig. 0Ja), which 
is in accordance with the prediction of Fig. 0Jb) . 

The discrepancy between the experimental data and 
the theoretical curve for weak lattices could be partly 
attributed to the small background of thermal atoms. 
These thermal atoms act as a reservoir as well, and cause 
decay through the Landau damping process (see, e.g. 
[9]). This results in somewhat more decay than theo- 
retically expected. For deep lattices, many-body effects 
cause depletion of the condensate on the order of a few 
percent, which could result in similar deviations from Bo- 
goliubov theory. In addition, small non-adiabatic effects 
in this regime contribute to the discrepancy between the 
actual experiment and the theoretical prediction. 

This is the first observation of finite-time (AZE-like) 
decay enhancement of a coherent, many-body excitation 
decaying into its natural environment. QZE is much 
harder to observe since (as explained above) it may only 
occur in weak lattices (where the GR decay is nonzero) 
and there decay dynamics is hard to access. The results 
presented here attest to the universality of Eq. (0| and 
its adequacy for multipartite systems. We have demon- 
strated that very limited knowledge of the environment 
response (its correlation time t c ), suffices for controllable 
modification of decay in complex macroscopic systems, 
if this control is faster than t c . Finally, this suggests 
the possibility of studying the quantum-classical transi- 
tion and decoherence dynamics using long-lived coherent 
excitations, and may be extended to other systems and 
their environments. 
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